Interaction of a magnetized pulsar wind with its surroundings 
MHD simulations of Pulsar Wind Nebulae 

E. van der Swaluw 
Dublin Institute for Advanced Studies, 5 Merrion Square, Dublin 2, Ireland 



O 
O 



Abstract. Magnetohydrodynamical simulations are presented of a magnetized pulsar wind interacting directly with 
the interstellar medium, or, in the case of a surrounding supernova remnant, with the associated freely expanding 
ejecta of the progenitor star. In both cases the simulations show that the pulsar wind nebula will be elongated due 
to the dynamical influence of the toroidal magnetic fields, which confirm predictions from a semi-analytical model 
presented by Begelman & Li. The simulations follow the expansion of the pulsar wind nebula when the latter is 
bounded by a strong shock and show that the expansion can be modeled with a standard power-law expansion 
rate. By performing different simulations with different magnetization parameters, I show that the latter weakly 
correlates with the elongation of the pulsar wind nebula. The results from the simulations are applied to determine 
the nature of the expansion rate of the pulsar wind nebula 3C58. It is shown that there is both observational and 
theoretical evidence which supports the scenario in which the pulsar wind nebula 3C58 has caught up with the 
reverse shock of the associated (but undetected) supernova remnant. 
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1. Introduction 

After the explosion of a massive star as a supernova, the 
progenitor star collapsed into a compact object, which can 
be a neutron star or a black hole. If the neutron star man- 
ifests itself as a rapidly rotating pulsar, the spin-down 
energy of the pulsar is converted into a relativistic pulsar 
wind (Rees & Gunn 1974), which blows a bubble or pulsar 
wind nebula (P WN) into the surrounding medium. The 
observational counterpart of a P WN is a plerion, which 
is characterized by filled-center diffuse emission (Weiler & 
Panagia 1978). For young plerionic supernova remnants, 
the central region contains the active young pulsar. 

The supernova explosion itself drives a blastwave into 
the interstellar medium (ISM), which results in an ex- 
panding shell of hot (shocked) material. The observational 
counterpart is a shell type supernova remnant, which 
is characterized by extended emission from the roughly 
spherical shell bounding the supernova remnant (SNR). 
Due to the deceleration of the expanding SNR by the 
surrounding ISM, a reverse shock forms (McKee 1974), 
which will propagate into the interior of the SNR when the 
amount of swept-up material becomes comparable with 
the ejected mass of the progenitor star. In the case of 
a SNR which contains both, an expanding shell of hot 
shocked material and a centered PWN (which is the obser- 
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vational counterpart of a composite remnant) , the interac- 
tion between both components becomes significant when 
the reverse shock hits the PWN (Reynolds & Chevalier 
1984; van der Swaluw et al. 2001; Blondin et al. 2001). 
After the passage of the reverse shock the PWN will make 
a transition from a supersonic to a subsonic expansion 
(van der Swaluw et al. 2001). The PWN will be deformed 
into a bow shock nebula when the motion of the pulsar 
becomes supersonic as the pulsar approaches the shell of 
the SNR (van der Swaluw et al. 2003). 

In this paper I will focus on young PWNe which are 
expanding supersonically, where the surrounding medium 
can be either the ISM or the freely expanding ejecta of 
the progenitor star (i.e. the stage of the PWN before it 
encounters the reverse shock) . For both scenarios there ex- 
ists observational evidence. In the case of SNR Gil. 2-0. 3 
a large PWN radius with respect to the shell of the SNR 
is detected at radio frequencies (Tarn et al. 2002), making 
the system a probable candidate for a PWN which has 
not encountered the reverse shock yet. In the case of the 
Crab Nebula there is indirect evidence for the presence 
of a shell containing shocked material and fast moving 
ejecta around the pulsar wind nebula (Sankrit & Hester 
1997). There are some counter-examples, naked plerions 
like G74.9+1.2 and G63. 7+1.1, where there is no evidence 
for any shell or fast-moving ejecta around the plerionic 
component (see e.g. Wallace et al. 1997a& b). This sug- 



gcsts that the pulsar wind m these systems is interacting 
directly with the surrounding ISM. 

A common feature of most plerions is a morphol- 
ogy with a significant elongation. This elongation is most 
strongly apparent at radio frequencies, where the complete 
bubble can be observed because of the long lifetime of ra- 
dio electrons with respect to the age of the PWN itself. 
Begelman & Li (1992) (hereafter B&L) explain the elonga- 
tion of the Crab Nebula by considering an axially symmet- 
ric model, where the elongation results from the dynamical 
influence of the toroidal magnetic fields in the pulsar wind 
bubble. Their model is a steady-state model, in contrast 
with recent time-dependent models of PWNe evolving in- 
side SNRs as discussed by van der Swaluw et al. (2001) 
and Blondin et al. (2001). However, these authors per- 
formed hydrodynamical simulations, and thus could not 
consider the dynamical influence of the magnetic fields on 
the PWN evolution. 

In this paper I consider the time-dependent evolution 
of a supersonically expanding PWN, driven by a mag- 
netized pulsar wind. The evolution of the PWN is mod- 
eled by performing magnetohydrodynamical (MHD) simu- 
lations using the Versatile Advection Code 1 (VAC). I con- 
firm the earlier analytical results from B&L, which state 
that the elongated morphology of the Crab Nebula (and 
other plerions) can be explained by the dynamical influ- 
ence of the toroidal magnetic fields in the pulsar wind 
bubble. 

I trace the position of the shock during the evolution 
of the PWN and show that although the magnetic fields 
are dynamically important, the expansion rate can still 
be modeled by a power-law expansion. I show that during 
the overall expansion the elongation will increase at a very 
slow rate, and it is at the very early stage of the PWN that 
its shape will be significantly elongated. Furthermore I will 
discuss the importance of the effects of the magnetization 
parameter <j, i.e. the standard ratio of the Poynting flux 
to the plasma kinetic energy flux, in the pulsar wind on 
the evolution of the PWN. Finally I will discuss the im- 
plications of the current simulations on the interpretation 
of the evolutionary status of the PWN 3C58. 

2. Evolution of a pulsar wind bubble 

2.1. The hydrodynamical limit 

I consider the evolution of a supersonically expanding 
PWN. This type of expansion occurs when the pulsar wind 
is directly interacting with the ISM or with the freely ex- 
panding ejecta of a young SNR. In this section I consider 
the hydrodynamical limit, i.e. no magnetic fields are in- 
volved, which enables one to use simple power-law expan- 
sion rates for the radius of the PWN (i? pW n oc <"). 

In the case of a pulsar wind with a constant luminosity 
Lq, interacting with the ISM with a uniform density po, 

1 See http://www.phys.uu.nl/~toth/ 



the eqation lor the radius ol the PWN is similar to the 
equation for stellar wind bubbles (Castor et al. 1975), i.e. 

/r \i/5 

R pvm = c(^\ t*l\ (1) 

where is of the order of unity. In the case of a PWN 
expanding in the freely expanding ejecta of a SNR, the 
equation for the radius of the PWN was derived by van 
der Swaluw et al. (2001): 

i? P wn = C 1 V t oc i 6 / 5 , (2) 

here Lq is the spindown luminosity, Eq is the explosion 
energy of the SNR, C is of the order of unity and Vq is 
defined by: 

En = fo , (3) 

here M c j denotes the mass of the ejecta. 

The ratio TZ between the radius of a PWN expanding 
into the ISM (equation 1) with respect to a PWN ex- 
panding into freely expanding ejecta (equation 2) is of the 
order of unity for supersonically expanding PWNe (using 
C=C=1): 

ft- 2.36 (-^-J FiooooC' 

here V\q ooo is expressed in units of 10 000 km/sec, E 5 i 
is the eplosion energy in units of 10 51 ergs, uq is the hy- 
drogen number density of the ISM in units of cm~ 3 and 
ikyr is the age expressed in units of 1 000 years. The sim- 
ulations indeed show an enhanced expansion rate for the 
pulsar wind expanding into the ISM for the considered 
ages (ikyr — 1). This can be explained by the fact that the 
pulsar wind embedded in the freely expanding ejecta will 
start to blow its bubble into the very dense medium of the 
ejected material from the supernova explosion. However, 
as the ejecta is freely expanding the density will decrease 
as a function of time (p e] cx t~ 3 ) which will accelerate the 
overall expansion of the pulsar wind bubble (e.g. Chevalier 
(1982); van der Swaluw et al. (2001)). This behaviour is 
reflected in the equations above by inserting po ~ t~ 3 
into equation (1) which yields the power-law behaviour 
of equation (2): the expansion of the pulsar wind bubble 
starts at a slow rate due to the dense medium but accel- 
erates to high expansion rates as the density of the ejecta 
is uniformly decreasing. 

I will use the MHD simulations to investigate the ex- 
pansion rate of a magnetized PWN, and compare the ex- 
pansion rate with the above power-law expansion rates. 

2.2. A magnetized quasi-static pulsar wind nebula 

A magnetized pulsar wind is driven by the spin-down lu- 
minosity of a rapidly rotating pulsar. The pulsar wind has 



a high Lorcntz factor (in the case ot the CJrab 1 w > 1U 
(Gallant et al. 2002)) and is terminated by a strong MHD 
shock, where the wind pressure and the confining ram 
pressure balance. At the position of the wind termination 
shock R ts one can write the luminosity of the pulsar wind 
as a combination of particle energy and magnetic energy 
(Kennel & Coroniti 1984): 



L = 47TTln w R^ s mc 3 (l + a) 



(4) 



where n w is the proper density in the wind, m the mean 
mass per particle and a is the ratio of magnetic energy to 
particle energy, i.e.: 
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a = 



47m w r w mc 2 ' 



(5) 



where B is the magnetic field component, which I will 
assume to be purely toroidal: 



B{r, z)cf>, 



(6) 



using an axially symmetric configuration (r, <j>, z). 

The shock heats the incoming flow and yields a high 
sound speed, c s <~ c/\/3 downstream of the MHD shock 
(Kennel & Coroniti 1984). 

B&L described the interior of a PWN by considering a 
steady-state solution (d/dt = 0) for the ideal MHD equa- 
tions in an axially symmetric configuration. Due to the 
high sound speed in the PWN, pressure imbalances are 
quickly smoothed out, therefore one can write the mo- 
mentum equation as: 
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From these two equations one can immediately observe 
that the total pressure along the symmetry axis (the z- 
axis) will be constant (equation 8), whereas the total pres- 
sure along the radial axis will decrease; the pressure gradi- 
ent in the radial direction balances the magnetic pinching 
force (the RHS of equation 9) , therefore the total pressure 
along the r-axis will decrease (equation 9). B&L show that 
the resulting difference along the PWN shock in pressure 
between the interior and the exterior of the bubble will 
result in an elongation of the PWN, as is observed in the 
Crab Nebula today. 





Simulation 1 


Simulation 2 


Pulsar wind luminosity L (erg/s) 


1.2 x 10 38 


1.2 x 10 38 


ISM mass density p (g/cm 3 ) 


i.o x irr 25 


1.0 x 10~ 23 


Terminal Velocity (cm/sec) 


3.0 x 10 9 


3.0 x 10 9 


Magnetization parameter a 


io- 3 


io- 3 


Adiabatic Index 


4/3 


4/3 


Explosion energy SNR E (erg) 


- 


10 51 


Mass ejecta M ej (M a ) 


- 


3 


Running Time (years) 


1,100 


1,400 


Number of grid colls (in r-dircction) 


100 


100 


Number of grid cells (in z-dircction) 


320 


320 


Grid size (pc) (in r-dircction) 


2.5 


2.5 


Grid size (pc) (in z-dircction) 


8.0 


8.0 



3. MHD simulations 

3.1. MHD calculations of a pulsar wind 

The difference between the model presented in this paper 
and the model of B&L is that I integrate the complete 
set of MHD equations in the non-relativistic limit with a 
continuous (thermal and magnetic) energy and mass in- 
put in a region centered around the position of the pulsar 
to model the pulsar wind. This yields the complete evo- 
lution of a PWN, for which one can trace the expansion, 
elongation rate and the internal structure, without being 
restricted to a quasi-static model as in B&L. 

I use the Versatile Advcction Code (VAC; Toth 1996) 
to integrate the equations of ideal MHD: 
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(10) 



0; (11) 



d fl 



(8) 



(9) h 



— + V x (V x B) = 0, 

here h is the enthalpy 

P 
P 



=(<»2) 



(13) 



p is the mass density, V is the velocity, B is the magnetic 
field strength, e is the internal energy density and P is the 
thermal pressure. 

I solve the MHD equations on a uniform grid in an ax- 
ially symmetric configuration (r, <fr, z) in 2.5D. In a 2.5D 
simulation all the flow variables T{ (where i=r,</>,z) are 
solved as a function of r and z, which includes the dynam- 
ical influence of the flow variables in the 4> direction, but 
the flow variables in this direction are constant at a given 
r and z, i.e. !F^{r, z) — T^r, z). I use a Total- Variation- 
Diminishing Lax-Friedrich scheme (Toth & Odstrcil 1996) 
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0.0 0.5 1.0 1.5 2.0 2.5 
Radius (in parsec) 

Fig. 1. Gray-scale representation of the total pressure dis- 
tribution of an elongated PWN. 



to solve the above MHD equations in conservative form. 
The boundary conditions are continuous everywhere ex- 
cept at the axis of symmetry where (a)symmetric bound- 
ary conditions were used for (momentum along the sym- 
metry axis) density, total energy and the momentum along 
the r-axis. 
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0.0 0.5 1.0 1.5 2.0 2.5 

Radius (in parsec) 

Fig. 2. Ratio of magnetic pressure P mag n, to the thermal 
pressure P t h throughout the PWN. A cut has been taken 
along the r-axis for different values of Z. 

3.2. Initial conditions of the pulsar wind 

I simulate a magnetized pulsar wind by depositing mass M 
and energy L continuously in a few grid cells concentrated 
around the position of the pulsar. The terminal velocity of 
the pulsar wind is determined from these two parameters, 
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0.0 0.5 1.0 1.5 2.0 2.5 
Radius (in parsec) 

Fig. 3. Gray-scale representation of the toroidal magnetic 
field distribution, where the labels in the plot are in units 
of 1(T 4 Gauss. 
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Fig. 4. Logarithmic gray-scale representation of the den- 
sity distribution, where the labels in the plot are in units 
of 1(T 25 g/cm 3 . 



i.e. Woo = y 2i/M pwn . In all simulations performed in this 
paper the terminal velocity has a value much larger then 
all the other velocities of interest. 

I want to include the dynamical effect of the magnetic 
fields in the pulsar wind bubble, therefore part of the to- 



tal luminosity L to t, is released as magnetic energy i magn ) 
where the magnetization parameter a, enters the equa- 
tions: 

imagn = ^tot = B 2 /2. (14) 




Fig. 5. Distance between the minimum (open triangles) 
and maximum (open squares) distance between the shock 
and the position of the pulsar (denoted in the text by R~ 
and R + respectively). The solid lines are curves which 
follow a i? pwn oc i 3 / 5 law. 

Since we are considering a PWN which is much larger than 
the size R\ c of the light-cylinder of the pulsar involved, 
we can safely neglect poloidal magnetic fields, therefore I 
will consider a purely toroidal magnetic field component. 
Simple estimates (c.f. Contopoulos & Kazanas, 2002) give 
a scaling law for the ratio of the toroidal B$, and poloidal 
Bp magnetic field strength, 

B^/B p ~ R/R lc » 1, (15) 

which makes the above assumption a valid one. 

Furthermore the current calculations do not solve the 
problem of transforming a high er-wind into a low er-wind 
(see e.g. Lyubarsky & Kirk 2001), the simulation already 
starts with a low value for a as is appropriate at the MHD 
wind termination shock (Kennel & Coroniti 1984). 

3.3. A PWN in an uniform ISM 

The simulation of a magnetized pulsar wind interacting 
with a homogeneous ISM has been performed with pa- 
rameters as denoted in table 1 (simulation 1). The MHD 
code evolves till a corresponding age of t = 1, 100 years. 
Below I will discuss the profiles of the different flow vari- 
ables of the system at the end of the simulation. 

Figure 1 shows the total pressure (= thermal pressure 
plus magnetic pressure) profile. One can observe the pres- 
ence of a pressure gradient along the r-axis, which balances 
the magnetic pinching force as discussed in the previous 
section. 



Fig. 6. Ratio of R + to R~ , i.e. the elongation of the pul- 
sar wind nebula due to the presence of toroidal magnetic 
fields. 

Figure 2 shows the ratio of the magnetic pressure and 
the thermal pressure where a cut has been taken along 
the r-axis at different positions on the z-axis. The model 
from B&L gives a ratio of magnetic to gas pressure which 
is monotonically increasing as a function of the cylindri- 
cal radius r and is independent of the z coordinate. The 
model presented here allows for local maxima, which is 
also reflected in the magnetic field distribution shown in 
figure 3. 

The density profile in figure 4 clearly shows the elon- 
gation of the PWN, where the ratio between the minimum 
and maximum distance between the shock bounding the 
PWN and the position of the pulsar as a function of time 
(denoted by R~ and R + respectively), equals ~ 1.5. This 
results confirm the elongation of a PWN by the dynamical 
influence of toroidal magnetic fields inside the pulsar wind 
bubble as predicted by the steady-state model of B&L. 

Figure 5 shows the evolution of the minimum and max- 
imum distance between the shock bounding the PWN and 
the position of the pulsar as a function of time (R~ and 
R + ). It is shown that the expansion rate for these shock 
positions still follow a simple power-law expansion rate, 
i.e. R + /~ oc t a with a = 3/5 as is appropriate for the 
expansion law of a purely hydrodynamical (stellar) wind 
interacting with a uniform ISM. This indicates that al- 
though the magnetic field influences the dynamics of the 
PWN, the evolution of a PWN driven by a magnetized 
pulsar wind can still be modeled by a standard power-law 
expansion rate. 

Figure 6 denotes the ratio between R + and R~ , from 
which one can observe that a significant part of the elon- 
gation of PWNe takes place in early stages (~ 100 years) 



ot its evolution. At later stages the elongation saturates at 
a very slow rate. In the last subsection I will investigate 
the dependency of the elongation on the magnetization 
parameter a. 

3.4. A PWN in freely expanding ejecta 

In order to investigate the structures observed in the sim- 
ulations presented above are dependent on the environ- 
ment in which the PWN expands, I performed a MHD 
simulation where a magnetized pulsar wind interacts with 
the freely expanding ejecta of a supernova remnant. The 
parameters of the simulation are denoted in table 1 (simu- 
lation 2). The freely expanding ejecta are initialised by im- 
pulsively releasing the mechanical energy and the ejected 
mass corresponding to a supernova explosion in the first 
few grid cells located around the position of the pulsar. 
Due to the larger scales of the expanding SNR with re- 
spect to the PWN, both the reverse shock and the forward 
shock will move off the grid. The freely expanding ejecta 
however will be represented correctly by using continuous 
boundary conditions, which means that the gradient is 
kept zero by copying the values of the flow variables from 
the edge of the grid into the ghost cells located around the 
grid. Using the above initialisation of the freely expanding 
ejecta of the SNR plus simultaneously simulating a pulsar 
wind as described above one can investigate the dynamics 
of the same pulsar wind as performed in simulation 1, but 
interacting with a different environment. 

At the end of the simulation, at a corresponding age of 
1, 400 years, the PWN has again been elongated along the 
symmetry axis. The overall structure of the different flow 
variables are similar to the ones discussed in the section 
above when the pulsar wind was interacting with a uni- 
form ISM. As an illustration figure 7 depicts the gray-scale 
representation of the magnetic field distribution which has 
the same overall structure as depicted in figure 3 for the 
first simulation. 

From the above results I conclude that a magnetized 
pulsar wind in either a uniform ISM or in freely expand- 
ing ejecta bounded by an expanding SNR, will yield an 
elongated PWN with a similar distribution of the several 
flow variables. This seems to imply that the shape and 
internal dynamics of a PWN is more determined by the 
central engine of the system then the surrounding medium 
in which it expands as long as the PWN is bounded by a 
shock. 

3.5. Elongation of pulsar wind nebulae 

The magnetization parameter determines the elongation 
rate of a supersonically expanding PWN. In order to quan- 
titatively determine the influence of this parameter I per- 
formed several simulations for a magnetized pulsar wind 
with parameters as in table 1 (simulation 1) but with dif- 
ferent values for the magnetization parameter. All values 
are in the range of earlier models of PWNc (Kcnnel& 




Radius (in parsec) 



Fig. 7. Gray-scale representation of the toroidal magnetic 
field distribution, where the labels in the plot are in units 
of 10~ 4 Gauss. 

Coroniti 1984; Emmering & Chevalier 1987; B& L). Figure 
8 summarises the results, which shows the elongation of 
the PWN as a function of time for the different MHD 
simulations. One can observe a weak dependency of the 
elongation on the magnetization parameter; an increase of 
roughly <~ 1.44 in elongation for an increase of 10 in the 
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Fig. 8. Elongation of a PWN driven by the same pulsar 
wind (parameters are denoted in table 1 (simulation 1)) 
but with different values for the magnetization parameter 
a. 



magnetization parameter. Furthermore, one can observe 
in all cases that the elongation saturates at late times in 
the evolution of the PWN. 



4. blvlK j(_bo, a crushed rWN ( 

A beautiful example of a plcrionic SNR which fits the mor- 
phology of the MHD results as discussed above is 3C58. 
The shape of this remnant shows a strong elongation at 
radio frequencies (see for example Bictcnholz et al. 2001). 
Recently the pulsar of this remnant has been detected, 
which places the active pulsar in the center of the rem- 
nant (Murray et al. 2002). This detection makes 3C58 
together with the Crab Nebula an excellent example of 
an elongated PWN blown into its surrounding medium by 
an active pulsar wind with the associated pulsar detected. 
I have shown in this paper that although the magnetic 
fields in the bubble are dynamically important and neces- 
sary to explain the elongation, the expansion of the PWN 
can still be modeled with a power-law expansion rate. This 
means that independent of the density of the surrounding 
medium, the expansion velocity of the PWN will roughly 
equal V^ xp = ai? pwn /t ag0 , where a = 6/5 when expansion 
takes place in the freely expanding ejecta of the progen- 
itor star and a = 3/5 when expansion is in an uniform 
ISM. Recent measurements of the radio expansion veloc- 
ity of 3C58 (Bietenholz et al. 2001) yield a value of ~ 900 
km/sec, which is in contrast with the observed values for 
-Rpwn = 4.5 pc and i age = 821 years (assuming an associa- 
tion with the historical SN event of 1181 and a power-law 
expansion rate), which would suggest an expansion veloc- 
ity of ~ 5000 km/sec. 

One way out of this discrepancy is to suggest that the 
remnant is much older. However this would suggest that 
the pulsar wind is interacting directly with the ISM, and 
there was no freely expanding envelope due to a very low- 
energy supernova event (Bietenholz et al. 2001). If this 
scenario represents the evolutionary status of 3C58, the 
shape of 3C58 corresponds exactly with the shape from 
the dynamical models presented in this paper. 

Another scenario for 3C58 was given by van dor 
Swaluw et al (2001), where it is suggested that the PWN 
of 3C58 is being compressed by the reverse shock. This 
scenario can explain a number of observational properties 
of 3C58 and would still associate the remnant with the 
presupernova progenitor of SN 1181: 

— The compression of the PWN by the reverse shock 
would decrease the expansion velocity, which is too low 
in the case of 3C58 if it was supersonically expanding 
(Bietenholz et al. 2001). 

— The compression of the PWN increases the magnetic 
field and particle energies, this can lead to an increase 
of the radio flux density as is observed in 3C58 in the 
period between 1967 and 1986 (Green 1987). 

— The compression of the PWN can bring down the fre- 
quency of the synchrotron cooling break (Gallant et al. 
2002), which is at a comparatively low frequency for 
3C58 (see e.g. Woltjer et al. 1997). 

— The passage of the reverse shock and the subsequent 
compression would also increase the pressure inside the 
PWN. This would decrease the radius of the wind ter- 
mination shock R ts - The ratio between the termination 



shock itts and the boundary oi the FW1N H pwn is in- 
deed remarkably small in 3C58, i.e. Rt s /R P wn ~ 0.01 
(Frail & Moffett 1993). 

If the above scenario is correct, the absence of any de- 
tected shell corresponding with the blastwave of the SNR 
remains a puzzle. 

5. Conclusions 

I have elaborated on the work of B&L by performing ax- 
ially symmetric 2.5D MHD simulations of a magnetized 
pulsar wind interacting with either a uniform interstellar 
medium or freely expanding ejecta. Both simulations were 
performed while the pulsar wind bubble is bounded by a 
strong shock. 

The simulations only included a toroidal magnetic field 
component, which was injected around the position of the 
pulsar. The simulations show an overall decline of the to- 
tal pressure in the radial direction in order to balance 
the magnetic pinching force due to the presence of these 
toroidal magnetic fields. The simulations show that the 
PWN is elongated as expected from the pressure differ- 
ences along the PWN shock. Furthermore the simulations 
show that significant elongation occurs at early stages of 
the expansion, while at later stages the elongation satu- 
rates. Finally by performing several simulations of a pulsar 
wind with different values for the magnetization parame- 
ter, I conclude that the elongation depends weakly on the 
value of the magnetization parameter. 

Although the simulations clearly show the dynamical 
importance of the magnetic fields, I have shown that one 
can still use a simple power-law expansion rate (i? pW n °c 
t a ) to model the expansion of the PWN. 

I have applied the implications of the current simula- 
tions on the naked plerion 3C58, which is an example of 
a PWN driven by a pulsar wind with a detected pulsar in 
the center. The morphology of this plerion fits perfectly 
well with the simulations performed for a magnetized pul- 
sar wind interacting with the ISM. However this would 
imply that the age of the remnant is much older and can 
not be associated with the historical supernova event in 
1181. Alternatively there is observational evidence that 
3C58 might be an example of a PWN being compressed 
by the reverse shock of the associated SNR. This would 
safeguard the association with the presupernova progeni- 
tor of SN 1181. Therefore in a future paper I will inves- 
tigate the physics of a magnetized pulsar wind embedded 
in an expanding SNR. In such a system the reverse shock 
will interact with an elongated PWN and would show the 
evolving plerionic morphology. 
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